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Theoretical studies of N 2 -broadened half-widths of H 2 0 lines involving high j states 
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(Received 20 September 2011; final version received 21 November 2011) 

Based on the properties of the energy levels and wave functions of H 2 0 states, one can categorize H 2 0 lines into 
individually defined groups such that within the same group, the energy levels and the wave functions associated 
with two paired lines have an identity property while those associated with different pairs have a similarity 
property. Meanwhile, by thoroughly analyzing processes used to calculate N 2 -broadened half-widths, it was 
found that the ‘Fourier series’ of (t\ j/Tj) and K , ( t ; y,T,), and a factor P 22 %(jf z fh x i) are the key items 

in the Robert-Bonamy formalism to distinguish contributions to ReS 2 (r c ) among different transitions of j/Tf<- jj 
T r However, these items are completely determined by the energy levels and the wave functions associated with 
their initial and final states and they must bear the latter’s features as well. Thus, it becomes obvious that for two 
paired lines in the same group, their calculated half-widths must be almost identical and the values associated 
with different pairs must vary smoothly as their values vary. Thus, the pair identity and the smooth variation 
rules are established within individual groups of lines. One can use these rules to screen half-width data listed in 
HITRAN and to improve the data accuracies. 

Keywords: energy levels and wave functions of H 2 0 states; Robert-Bonamy theory; theoretically calculated 
half-widths and shifts; rules governing spectroscopic parameters of H 2 0 lines; HITRAN 


1. Introduction 

The modelling of the atmosphere from satellite-based, 
balloon-based, and Earth-based instruments requires 
an accurate spectroscopic database such as HITRAN 

[1.2] . This widely used database has spectroscopic 
parameters for the most important molecules in bands 
from the microwave to the ultraviolet spectral regions. 
It is obvious that the accuracy of these data is essential 
for users’ applications. This is especially true in 
accurate atmospheric retrievals involving the very 
important water vapor molecule. In order to meet the 
accuracy requirement for the H 2 0 molecule, the H 2 0 
database in HITRAN has been updated several times 

[1.2] and a new HITRAN 2012 version will appear 
soon. In the recent updating process, the majority of 
the pressure broadened half-widths, the temperature 
exponents, and the pressure-induced shifts come from 
theoretically calculated values using the Robert- 
Bonamy (RB) formalism [3]. It is worth mentioning 
that the RB formalism was developed more than three 
decades ago and it contains several basic assumptions 
whose applicability was not thoroughly justified. 
Mainly due to difficulties in developing better 


formalisms to carry out practical calculations, there 
are no viable alternatives available and one has to rely 
on the RB formalism at present. Besides, for lines with 
small half-width values, contributions to their calcu- 
lated half-widths result mainly from nearly head-on 
collisions in which two interacting molecules can reach 
closer distances [4]. Unfortunately, the largest uncer- 
tainties in modeling interaction potentials and collision 
trajectories usually occur in the short distance regions 
and they could lead to large errors in calculating 
contributions from nearly head-on collisions. Thus, 
calculated half-width values for lines with small half- 
width values are less reliable [4]. Meanwhile it turns out 
that, in general, lines with small half-widths are those 
involving high j states. Thus, uncertainties associated 
with calculated values would be worse than the 
requirements needed for users’ atmospheric applica- 
tions, especially for lines with high / values. 

On the other hand, some of the data for the water 
vapor lines are experimental and may come from 
different labs. Given the fact that lines with high j 
values are usually weak, to perform measurements for 
them becomes more difficult and measured results 
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could contain larger uncertainties. In addition, there 
are a huge number lines especially at high tempera- 
tures, and the majority of them are associated with 
high j states, and thus to provide these data from 
experimental measurements is not realistic. In sum- 
mary, the current uncertainties for the half-widths in 
the H 2 0 database are far beyond the desired accuracy 
and it seems that this situation could last for a while. 

In order to find realistic ways to improve the 
accuracies of half-widths of H 2 0 lines, we began this 
study by wondering whether links between lines of 
interest and their half-widths do exist or not. If the 
answer is yes, we hoped that some rules representing 
these links can be discovered. Of course, the success of 
this idea relies not only on the nature of pressure 
broadened half-widths, but also the answer to the 
fundamental question about why and how calculated 
N 2 -broadened half-widths vary for the lines of interest. 
In fact, we have been interested in this fundamental 
question for several years and have made many efforts 
to search for the answer. Very recently, however, we 
have found a method and in the following we outline 
our route leading to the answer yes. Readers who are 
not familiar with the theory may skip over the next three 
paragraphs. 

It is well known that calculated half-widths depend 
on the lines of interest and also on the environment in 
which the absorber molecule H 2 0 is immersed. 
In order to find the answer, our approach is to 
separate these two dependences as completely as 
possible. In theoretical calculations, the former is 
described by the initial and final states of H 2 0 and 
the latter is explicitly represented by the interaction 
potential and trajectory models. Within the RB 
formalism, our investigation targets are contributions 
to the term of Re5 2 (r c .) from individual correlation 
functions because Re,S' 2 (r r ) plays a crucial role in 
determining calculated half-widths. Along this way, we 
have found that for each of the correlations, both its 
contributions to ReS 2 , oute r, /(a) and R eS 2 , outer, f(r c ), 
two of the components of Rc,S' 2 (/y), can be expressed as 
integrations over time whose integrands mainly consist 
of two functions. One is the correlation function 
FuK^Liit) that depends on the potential and trajec- 
tory models and is common for all lines. The other is 
the two functions K ,(t\ jx) which depend on the 

initial or the final states oif the lines of interest but are 
independent of the environment. Thus, for an individ- 
ual line, two If^j v (t;jr) associated with its initial 
and final states contains all information necessary to 
determine its contributions to Re.S' 2i outer ,{r c ) and 
R eS 2 , outer, f( r c)- Therefore, it is these two functions 
^LiK\K / (0 that distinguish the different magnitudes of 
R eS 2 , outer, i(r c ) and R eS 2 , outer, f(r c ) for individual lines. 


Meanwhile, similar analysis can be applied to the term 
of S 2 ,middie( r c), the remaining component of R eS 2 (r c ) 
also. In this case, it is a factor P that 
depends on the lines of interest but is independent of 
the environment. Then, we know that in order to find 
the answers we need to focus our attention on the two 
functions (t;jr) and the factor P Lx K x K^jfifji*i)- 

It turns out that expressions for v (/; /tt) are 

given in terms of ‘Fourier series’ and the number of 
their major components is very limited. This implies 
that each \ VL K K ,{t\jx) can be well represented by a 
spectrum containing only several components. As a 
result, to analyze profiles of ^ (t;jr) among many 
different states becomes easier, hy simply plotting 
distributions of the spectra associated with different 
H 2 0 states, we have found that there are identities of 
their distributions between certain pairs of states and 
there are similarities within certain groups of states 
also. With respect to values of P L,K,K\(jf'Cfji'Ci), similar 
features have also been found. In addition, the events 
associated with R^ A| ) A|A '(C// T /) and ^ (fi ./A,) and 
the events with Pr,,K, k\ \jf Tfjiti) happen simultaneously 
among the same lines. 

Thus, we are sure that the similar features must 
exist for the calculated half-widths among certain lines. 
Armed with the knowledge where the identity and 
the similarity happen for K ,(t; jr) and 

Pl, k, k\ UfTjj,?,), we know how to categorize lines into 
different groups. Then by checking values of the 
calculated half-width for lines in the same group, our 
expectation is well verified. At this stage, we are only 
one step away from finding the answer, namely, we 
have to explain why these features exist within the 
groups. Based on quantum mechanics, the energy 
levels and the wave functions of the H 2 0 states contain 
all information about the states. Therefore, the latter is 
a right place to search the answer. By investigating the 
properties of the energy levels and the wave functions 
of states, their identities and similarities clearly emerge 
within the groups. Thus, we have got the answers. With 
respect to the first question why calculated half-widths 
vary with the lines, it results from individual energy 
levels and wave functions associated with their initial 
and final states. Concerning the second question how 
calculated results vary with the lines, we know the 
answer within individually defined groups. It is the 
identities and the similarities of their energy levels and 
their wave functions that govern their variations. 

After realizing the above, an immediate thought 
was that there must be similar rules valid for other 
spectroscopic parameters. In fact, we consider a system 
consisting of an absorber water molecule immersed in 
bath molecules while the absorber interacts with 
electromagnetic fields as a black box. The inputs and 
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the outputs of this black box are the transition lines of 
H 2 0 and their spectroscopic parameters. The main 
assumptions introduced here are: (1) the outputs 
depend on the inputs; (2) identical inputs should 
yield identical outputs; and (3) similar inputs should 
yield similar outputs. Then, by fully exploiting the 
identity and similarity properties of the energy levels 
and wave functions, without knowing what really 
happens inside the black box, we can expect that the 
identity and similarity features exist among the spec- 
troscopic parameters. We have found that measure- 
ments from different groups do demonstrate the above 
claim. In addition, several important applications have 
been proposed. The study based on the black box 
theory was summarized in a recent paper [5], 

For the theoretically calculated N 2 -broadened half- 
widths, we know completely what happens inside of 
the system. Instead of relying on the black box, we can 
establish these rules from rigorous theoretical analyses. 
In addition by developing tools used in the analysis, 
new features of the line-shape theory can be exhibited. 
Therefore, the present work is supplementary to recent 
work [5], In the present study, we do not follow our 
discovery route exactly. Because properties of the 
energy levels and wave functions are fundamental, we 
first present our discussion of this subject in Section 2. 
With this as a guide, we then present our analysis, step 
by step, in the following sections. Finally, we present a 
short discussion and conclusion section. 


2. Properties of energy levels and 
wave functions of H 2 0 

It is well known that when one studies energy levels and 
wave functions for an asymmetric top molecule such as 
FI 2 0, one must identify its three principal axes a, h, c 
with a body-fixed frame [6-9]. There are six different 
choices to make the identification and they are 
designated: I R, I L, II R, II L, III R, and III L. Of 
course, in describing energy levels and wave functions, 
all these representations are equivalent. The energy 
levels are independent of the choice of representation. 
Flowever, the wave functions would be different, but 
wave functions derived from different representation 
are exchangeable. One can change them by using 
rotational transformations representing proper rotation 
operators under which one body-fixed frame rotates 
into another. In practice, one can choose the most 
suitable representation according to one’s preference. 

For years there have existed subroutines available to 
enable one to obtain energy levels and wave functions 
for important rotational states. Thus, in pursuing their 
spectroscopic studies, most people adopt energy levels 


and wave functions derived from these subroutines and 
they usually do not pay too much attention to these 
items themselves. However, it turns out that exploiting 
the properties of energy levels and wave functions 
proves to be a great help. In order to support this claim, 
we exhibit some properties of the energy levels and wave 
functions of H 2 0 in the ground state first. 

Before continuing our discussions, we note that 
usually states of the H 2 0 molecule are labelled by three 
quantum numbers j, k a , and k c . In the present study, 
sometimes we may use j and r ( =k a — k c ) to label a 
state because we want to use short labels, especially in 
figures. There is a simple way to recover the corre- 
sponding k a and k c values from j and r. For the case 
where j+r— 2n, k a — n and k c —j—n. Meanwhile, for 
the case where j + z — 2n + 1, k a — n + 1 and k c — j — n. 


2.1. Properties of energy levels 

With respect to the energy levels, we present rotational 
energy levels of H 2 0 with j = 1 1 to 20 in the vibrational 
ground state provided by Barber et al. [10] versus j in 
Figure 1. In this plot energy levels of states with 
j + k a — k c — even are plotted by symbols x and those 
with j + k a — k c — odd are given by symbols A, respec- 
tively. In addition, for the former, their values of 
k a — k c are presented on the right side of x and for the 
latter, the values are on the left side of A. As shown in 
the figure, there are identities of the energy levels 
between pairs of the states with the same j. Thus, 
within certain small tolerances, the pair identity of the 
energy levels holds well. The higher the j is, the firmer 
the identity holds. With Figure 1, one can easily 
identify those paired states. They are pairs of two states 
with the same j and k c (—j — n), but their k a — n or n + 1 
where n — 0, 1, . . . or two states with the same j and k a 
( —j — n ), but k c — n or n+ 1. We list some of them 
in braces: {jij-ujij-i}, {jjj-Ufj-i}, ■ ■ ■ , 

{ jj-2,3,jj-2,2},{jj- i, 2 j)- 1 , i } : and {jj,ujj, o}- Besides this 
pair identity feature, one can conclude that for each of 
these pairs, their energy levels vary smoothly as j varies 
and these variation patterns are well organized along 
the j axis. 


2.2. Properties of wave functions 

With respect to wave functions of the H 2 0 states, they 
are given in terms of expansion coefficients U J kr over 
the symmetric top wave functions jkm) with k — —j, 
~j+ 1 j h 

\jxm) = ^2 U iz \jkm). (1) 

k 
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Figure 1. A plot to show energy levels of H 2 0 states with j= 11 to 20 in the vibrational ground state [10], For states with 
j + k a — k c = even, their energy levels are plotted by x and their values of k a — k c are presented on the right side of the symbols. 
Meanwhile, for states with j +k a — k c = odd, their energy levels are plotted by A and values of k a — k c are on the left side of the 
symbols. 


Table 1. Symmetry classification by the parities of k a and k c . 



1 R 

I L 

II R 

11 L 

III R 

III L 

J even 

J odd 

*^even 

J odd 

J even 

Jodd 

Jeven 

Jodd 

J even 

Jodd 

J even 

J odd 

E + 

ee 

eo 

ee 

eo 

ee 

OO 

ee 

OO 

ee 

oe 

ee 

oe 

E" 

eo 

ee 

eo 

ee 

00 

ee 

OO 

ee 

oe 

ee 

oe 

ee 

o + 

00 

oe 

oe 

00 

oe 

eo 

eo 

oe 

eo 

OO 

OO 

eo 

cr 

oe 

00 

00 

oe 

eo 

oe 

oe 

eo 

OO 

eo 

eo 

OO 


With the selected representation, states assigned either 
by the sub-block E + or by E“ have non-zero values of 
U kz only for k is even and states associated with 0 + or 
O - have non-zero values only for k is odd. In addition, 
depending on the superscripts + and — in the sub- 
block symbols assigned, U kz is either an even function 
of k or an odd function of k. In the present study, we 
follow a convention adopted by Zare [9] to define the 
Wang functions and the symmetric top wave functions. 
With this convention, symmetry classification by the 
parities of k a and k c in the six representations is given 
by Table 1. 

After outlining the basics about how to describe the 
wave functions, we exhibit their properties. First, we 
consider three sets of pairs of states {jj,o,jj, i}, 
0 }_ u ,y}_ i j2 }, and {;}- 2 , 2 j}- 2 , 3 }and present their coef- 
ficients U kz starting from certain boundary values j bd 
(i.e., 3, 5, and 7, respectively) to j— 26 in corresponding 


Figures 2(a)-(c). For these states whose k a are the 
maximum or close to the maximum, the I R represen- 
tation is the best choice. 

As shown in the figure, there are two striking 
features of these coefficients. First of all, coefficients of 
the paired states such as j J0 and 7 }^ with the same / are 
almost identical in the positive k axis and one becomes 
the other’s mirror image in the negative k axis. In other 
words, their absolute values are always the same. It is 
worth mentioning that their opposite signs in the 
negative k axis are necessary to guarantee the orthog- 
onality between their wave functions. Neglecting the 
sign difference in the negative k axis, we call this 
feature the pair identity of the wave functions. We note 
that the higher the j, the better the pair identity. 
Because the pair identity breaks down for j < / bd , the 
plotting in Figure 2 starts from j — j b d . Secondly, by 
comparing the coefficients associated with different 
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Figure 2. A plot to show properties of H 2 0 wave functions for three sets of pairs of states {y’/,o J/,i }, {.//— i , i , J/— [ , 2 } , and {jj-2,2, 
jj- 2 , 1 } for j~Jbi, ■ ■ ■ ,26. Their coefficients U J kr with k ranging from —j to j in the 1 R representation are represented by 
symbols x and A, respectively. Different colors are used to distinguish different j values and symbols associated with the same 
states are connected by dotted or solid lines. 


j values, it is obvious that for members of the pairs 
such as jj_ 0 with j— 3,4, , their patterns are very 
similar. As j varies, without any distortions, place- 
ments of the patterns shift horizontally along both the 
positive and negative k axes. We call this the pattern 
similarity of the wave functions. It is worthwhile 
emphasizing that these pattern similarities are valid 
only among members of these pairs. For example, they 
are valid within a set of the states jj 0 and within a set 
of the states y) 1 , and so on. Finally, by comparing 
Figures 2(a)-(c), values of the boundary y'bd increase 
from 3 to 5 to 7, respectively. This implies that as k a 
decreases from the maximum, the pair identity and the 
smooth variation break down at higher and higher 
j values. 


Now, we consider states whose k c values are the 
maximum or close to it. In these cases, the III R 
representation is more suitable to describe the proper- 
ties of the wave functions. We present coefficients U J kz 
for three sets of pairs of states {j 0 j,jij}, {jij-iJij-i}, 
and {J2J-2J3J-2} starting from certain boundary 
values y' b d (be., 7, 10, and 13, respectively) to y = 26 in 
Figure 3. As shown in the figure, the pair identities and 
the pattern similarities appearing in Figure 2 for the 
previous three sets of pairs clearly hold there also, but 
it is noticeable that their validity boundaries become 
higher than those associated with the previous 
three sets. 

In practice, for sets of paired states, one can find 
boundary estimations from plots to show their wave 
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k 


Figure 3. The same as Figure 2 except that the three sets of pairs of states are { 70 , 7 , 71 , 7 }, {jij-iJij-i}* and { 72 , 7 - 2 , 73 , 7 - 2 } an d 
their wave functions are derived in the III R representation. 


Table 2. Boundaries for different sets of paired H 2 0 states. 


Set 

Jj\ o 

7)- i,i 

jj- 2.2 

jj- 3,3 

jj- 4,4 

jj- 5,5 

jj— 6,6 

74,7-4 

73J-3 

72,7-2 

71,7-1 

70,7 

Jj, 1 

Jj- 1,2 

jj- 2,3 

jj- 3,4 

jj- 4,5 

jj- 5,6 

jj- 6,7 

75,y-4 

74,7-3 

73,7-2 

72,7-1 

71.7 

jbd 

3 

5 

7 

9 

10 

12 

14 

19 

16 

13 

10 

7 

E (%) 

0.63 

0.66 

0.49 

0.32 

1.08 

0.60 

0.32 

1.28 

0.55 

1.13 

1.00 

0.60 


functions such as Figures 2 and 3. Alternatively, for a 
set of paired states j,i\ and y,r 2 , by defining the relative 
difference of their wave functions e as 


s — 




one can calculate how e varies with j and determine a 
boundary for this set. We note that e is a numerical 


measure of the identity of the paired wave functions. 
We provide our suggested j bd and their corresponding 
values of s in Table 2. As shown in the table, we have 
chosen s about 1% to determine y' bd and the latter 
match those estimated from the figures well. 

In summary, with respect to the FFO energy levels, 
there are pair identities and well organized variations 
exhibited within each of the sets of paired states whose 
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j values are high. Meanwhile, for the wave functions, 
there are pair identities and pattern similarities within 
each of the same sets defined in analysing the energy 
levels. In general, sets with small differences between j 
and \k a — k c \ have lower boundaries / bd and those with 
large differences have higher ones. On the other hand, 
except for the pair identities, other properties of the 
energy levels and wave functions have individual 
characters. For example, similar patterns of the wave 
functions appearing in one set differ significantly from 
those in other sets. Furthermore, as explained above, 
these properties disappear for states whose j values are 
below certain boundaries and the latter vary with 
the sets. 

There are useful explanations about symmetry 
assignments of H 2 0 wave functions for paired states. 
Besides, there are quantities associated with the 
energy levels and the wave functions which are 
important in the present study because they determine 
couplings between H 2 0 states. Readers can find 
discussions of these in sections A-l and A-2 of the 
Appendix. 

Finally, we would like briefly to point out that it is 
not difficult to explain why the paired states with high j 
values have almost the same energy levels and the wave 
functions. It is well known that the Hamiltonian of the 
H 2 0 molecule associated with the quantum number j 
can be expressed by a sub-block form consisting of 
four sub-matrices, E + , E“, 0 + , and O - . As an 
example, we assume y = even and consider the pairs 
of jjfi and jj \. According to Table A-l in the Appendix 
and Figure 1, is an eigenvector of E + with the 
largest eigenvalue and j- h \ is an eigenvector of E _ with 
the largest eigenvalue. Both of these E + and E“ are 
tridiagonal and they are almost identical except for 
adding additional first low and first column in E + . In 
addition, in comparison with their diagonal elements, 
their off-diagonal elements are small. The higher the j 
is, the more identical are E + and E“. It turns out that 
their largest eigenvectors and eigenvalues are their ones 
least affected by their differences. As a result, they bear 
the identity feature the most. For other cases, similar 
arguments are also applicable, but we do not pursue 
this here. 


3. General formalism in calculating half-widths for 
H 2 0 lines 

After discussing the properties of the energy levels, 
wave functions, we briefly outline the formalism used 
to calculate the half-widths and shifts of H 2 0 lines. 

The main computational task for calculating the 
Lorentzian half-widths is the evaluations of matrix 


elements appearing in the perturbation expansion of 
the S matrix (—Si • Sp, where Si and S F are scattering 
matrices in Hilbert space). Usually, in practice, these 
evaluations are limited to the second-order of the 
expansion. Within the RB formalism that has been 
widely used for calculating Lorentzian spectral line 
half-widths and shifts for decades [11-14], the original 
expression for the half-width is given by 

n h f + °° f + °° 

Yrb — z — / vf(v) dv Ijthdh 

C Jo Jo 

X ( 1 - costs) ( h ) + ImS 2 (h)]e- Re52(i,) ) 7 . , (3) 

where n h is the number density of the bath molecule, /(v) 
is the Maxwell-Boltzmann distribution function, Si and 
S 2 are matrix elements in the Liouville space associated 
with the first- and second-orders of the perturbation 
expansion of the Liouville operator S, and { ) J2 means 
an average over the quantum number j 2 of the N 2 
molecule. However, we have found a subtle error in the 
RB formalism [15]. After remedying this derivation 
error, the correct expression for the half-width in the 
‘modified’ RB formalism (MRB) becomes [15] 

n b r+°° r+°° 

Ymrb = ~ / v/'(v)dv / 2nb d h 

2 TCC J o Jo 

X { 1 - cos[(Si(b))j 2 + Im(S 2 (h)> A ]e- Re < S2 ®>A }. 

(4) 

We note that the essential difference between these two 
expressions is that in Equation (3) the summation over 
j 2 is outside of the cumulant expansion while in 
contrast, in Equation (4) it is inside. It is worthwhile 
emphasizing that because the bath average is carried 
out in the line space, an average over j 2 has been 
included. Thus, by understanding that Si(b) and S 2 (b) 
are the bath averages associated with the first- and 
second-order expansions of the Liouville operator, one 
can supress all of these ( ) /2 in Equation (4). For 
simplicity, we will omit the subscript MRB of y. 

In the present study, we consider lines in the H 2 0 
pure rotational band. Then, by ignoring negligible 
contributions from Im.S' 2 (b), making an approximation 
to replace the integration over the velocity by the 
averaged velocity v(= y/SkT/jrnt) and changing the 
integration variable impact parameter b and the lower 
limit 0 to the distance of closest approach r c along 
trajectories and r cmm , the minimum value of r c corre- 
sponding to strictly head-on collisions, respectively, one 
can obtain a simplified expression for the half-width 
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The value of r, min can be determined by an energy 
conservation equation 

2 Fisovc, min) in rn\ 

=2 1 = 0. (6) 

mxr 

Usually, people prefer to represent S 2 by three 
components labelled by S 2 , outer, i, S 2 , ou terj, and 
middle, respectively. In the present study, we follow 
the same custom. As an example, an expression for 
S 2 , outer, i where the subscript i is a simple notation for 
the initial state /,■ r ; is given by [4] 



x ^■x i w i j 2 m 2 \V(R(t))\j.x , i m , i j 2 m'^j 

x r / w /./2 w 2 I V(R(t ')) I jirmjjm^j. (7) 

In the above equation, pj 2 is the density matrix of 
the bath molecule, coj^j'r' is given by equation (A- 1 ), 
coj y' = [rt b \j 2 ) - d b \ j 2 )]/ti , and \jiXimij 2 m 2 ) = 

7,T/w,-) <g> \ j 2 m 2 ) where | /Vt//77,) and \j 2 m 2 ) are the 
basis of Hilbert space for the H 2 0 and N 2 molecules, 
respectively. 

In order to overcome convergence problems exist- 
ing in usual line-shape theories involving the site-site 
potential models, we have developed a new formalism 
based on the coordinate representation. We will not 
provide a detailed development of the new formalism, 
but interested readers can find its derivation in our 
previous work [4,16], Here, we only present some of 
the important formulas. By introducing the correlation 
functions F j A K t K\L 2 {t) and their Fourier transforms 
H L x K i K , l L 2 {oi) defined by 

— —f= / e m 'FL l K l K , l L 2 (t)dt. ( 8 ) 
V 2.7T J —oo 

the expression for the real part of S 2 oute , 2 can be 
obtained as 



x DiJiXijft, L\ K\ )D(j i Tij’ i T; ; L\ K\ ) 

x^(2 /2 +l)(2 / -' + l)p /2 

hi'i 

x C 2 (j 2 j 2 L 2 , 000) 

X H LlKlKi L 2 ((Oj i T i j' i r' t + w jzi 2 )- (9) 


We note that because the correlation functions 
FhKiiCLi.it) are associated with a specified trajectory, 
ReS 2j ou ter, i is a function of r c . For clarity, we have 
explicitly added the argument r c for Re.S' 2 rm/cr , ,■ in 
Equation (9). An expression for Re.S’ 2 oH ; cr j(/' c ) is the 
same as Re5 2j outer , ,(r c ) given in Equation (9) except for 
a replacement of j h r ( - by jf, x f . Meanwhile, an 
expression for S 2 , middle (r c ) is given by 

S 2 ,middie(r c ) = J In Y, {(-l) L| (2y/+ l)(2j/+l) 

L\ V, K\ L 2 

x mjujjijfi | 7 -i) 

X DijiXijitr, L i K\ )D(j f tfjf r f ; Li^)} 

xY(2j2 + \m+\)p j7 

jU'2 

x C 1 {j 2 j' 1 L 2 , 000)H LlKlK ' iL2 (Mj 2 jF). 

( 10 ) 

We note that with our new formalism, the main 
tasks to calculate N 2 -broadened half-widths for H 2 0 
lines are evaluations of several dozens of the correla- 
tions F LlKlK ' iL2 (t) labeled by a set of numbers consisting 
of one tensor rank L x with two subsidiary indices Aj, 
K\ related to H 2 0 and another tensor rank L 2 for N 2 
[4]. Because N 2 is a diatomic molecule, L 2 must be 
even. Thus, the number of sets is determined by the 
upper limits of L\ and L 2 . If one chooses the II R 
representation to develop the H 2 0 wave functions 
where the two H atoms are symmetrically located in 
the molecular-fixed frame, values of K\ and K\ must 
also be even. Due to symmetries, some of the corre- 
lation functions are identical. For examples, the four 
correlations labeled by (2 2 2 0), (2 2—2 0), (2—2 2 0) and 
(2—2—2 0) are identical. One can conveniently express 
this symmetry by F Li \k,\\k\\l 2 (0- In addition, there is 
another exchange symmetry between Aj and K\ in the 
correlation functions. Thus, the number of the corre- 
lations required to be evaluated can be reduced 
significantly. 

Finally, as shown in Equations (9) and (10), the 
expressions for R eS 2 ^ outerJ (r c ), R eS 2 ^ outer j(r c ), and 
S 2 , middieO'c) ar e given by summations of contributions 
from individual correlations. As a result, among all 
correlations how important the individual correlation 
is depends on how large its magnitude is. For a 
specified correlation, the asymptotic behaviour of its 
magnitude as r f ->oo is well defined. In general, 
(1000), (1002), (2|2||2|0), and (2|2||2|2) have large 
magnitudes and at r c — > 00 they vary as r'" with 
m — — 14, —8, —16, and —10, respectively. In the 
present study, we will focus our attention on these 
correlations. 
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4. Theoretical tools in analyzing calculated 
half-widths 

4.1. Harmonic expansions appearing in 

$2, outer, i(^c) Olid Sj outer, i( C t ) 

In the process of deriving Equation (9) for 
RcS 2 , outer, life) from Equation (7), there is an interme- 
diate expression for S 2> outer, ti r c) such that 


*~* 2 , outer, iv'c) 

poo 

= E / 

l,k i k\l 2 Jo 

( 11 ) 

where two functions are defined by 

- E (2 /" + 1W/T.7#; L ^) 

j'i z 'i 

xDUmj^-Li^e 10 ^, (12) 

and 

<V) = E^+ 1 )^2+ 1 H 


x C 2 (j£L 2 ,m)e m »l‘, (13) 


respectively. It is worth mentioning that the correlation 
functions F^KiKLiif) depend on the potential and 
trajectory models. In contrast, the functions 
k’ (0 j‘ T d and IE^(t) are independent of these. 
On the other hand, F^k^LiO) and are 

common for all lines. Meanwhile, the functions 
K'( l > 7 T ) associated with j, and r,- and their 
partners associated with jf and t y are only factors 
depending on the line of interest. They contain all 
information necessary to distinguish different amounts 
of contributions to the half-width of individual lines 
from the correlations to be considered. 

It turns out that in Equations (12) and (13), both 
and <V) are given in terms of 
harmonic expansions. Although these two functions 
are not periodic, their forms are very similar to the 
usual Fourier series. Thus, we borrow the concept of 
the Fourier series here. For later convenience, we 
rewrite Equation (12) as 


As usual, the ‘Fourier series’ of w K ’ (0 jAt) are 
given by a set of components consisting of two values 
of co a {jjti) and A a (L\K\E x \jjXj). The total number of 
components equals how many choices of a, and in 
general, there are several dozens of components. 
Fortunately, as shown later, there are only a few 
strong ones because many components are very weak. 
With Equations (A- 2) and (15), one can show that 

E A a (LiKiK l l -JiTi) = (2L| + 1)$**-. (16) 

O' 

The summation of A a {L\ K\ K \ ; y,r,) over all a equals to 
2L| + 1 if ATj =K\ or equals to 0 if K x / K\ . We note 
that the above discussion is also applicable for 
analysing ReS 2i outer, where the ‘Fourier series’ of 
^l,k, K' (Pjf x f) are introduced. 

Similarly, we can rewrite Equation (13) as 

IF®(0 = E s ^ 2 K^ (17) 

p 

where the summation index represents a selected pair 
of j2 and /j, a>p stands for Wj y -', and Bp(L 2 ) is given by 

Bp(Fi) = (2 72 + 1)(2 /' + \)p h C\j 2 j'2 L 2 , 000). (18) 

Meanwhile, it is easy to show that the total intensity of 
the ‘Fourier series’ of W^(t) is equal to 2L 2 + 1. In 
contrast with A a (L\K\ K \ ; y,r,) which are independent 
of the temperature, the ‘Fourier series’ of W)y(f) 
depends slightly on the temperature. 

Furthermore, by combining Equation (14) and 
Equation (17), we have 

<W^ r M A ,V) 

= EE^ a (L\K\ K \ ; jiTj)Bp(L2)e^ tt>a( j‘ T d +co d t 

a p 

= E£k(^i^2;.// xd^ Um)t , (19) 

Y 

where y is a simple notation for a set of j', r', j 2 , j 2 , 

C)y( j jti) — Olj.z.j’C T , aild Fy(l.\ K\ A| l.l'.j-Ti) — 

A a {L\K\K x \ jjXj)Bp(L 2 ). After having the expressions 
for E y (L\ K\ K\ LyJtXi) and co y ( jiXi) available and 
knowing all accessible choices of y, one can simply 
rewrite ReS 2 , outer, ,{r c ) as 


- E Aa(L ' K ' K'djAde mrAj,Ti> ‘, (14) 

a 

where the summation index a represents a selected 
combination of /■ and x- in Equation (12), co a (jiXi) 
stands for coj^pg, and A a (L\ K\ A") ; y,T,) is defined by 

A a (L\ K] K\ ;jiXi) 

- (2 % + VDUiXijft, LiKODUtXijft, FK\ ). (15) 


R eS^oVer/Tc) = J* E 

" - LtKtK^Lz y 

X H LiKi K \ Li [0>y (. jiTi)\ ■ (20) 

This expression clearly shows for an H 2 0 line of 
interest, how contributions to the ReS 2 outer, t( r c) term 
from individual correlations are calculated. First of all, 
one determines a set of pairs E Y (L\K\K { L 2 \jiXj) and 
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(Oy(jiTi) with all accessible choices of the y. We note 
that in contrast with the correlations that are common 
for all H 2 0 lines and depend on the interaction and 
trajectory models, this set of pairs depends on the line 
of interest but is independent of the interaction and 
trajectory models. In the next step, for an individual 
correlation, one picks up values of its Fourier trans- 
form k\ /. 2 (&>) at cd — coy(jiTi), multiplies the latter 
by the corresponding values of E y (L\ K\ K\ L 2 \j{ii), and 
adds up all the results to obtain the total contributions 
from this correlation. One then repeats the above step 
for each of the correlations, and finally one adds up 
contributions from all correlations to obtain the 
final result. 


4.2. General properties of the ‘Fourier series’ of 

We outline some properties of K f t\ jr) that are 

very helpful in calculating contributions to 
R es 2 , outer,, (rj and ReS 2 , outer, f(r c ) from the correla- 
tions. For simplicity, as long as their arguments are not 
involved in discussions, is simplified to 

W ^ K K ' ■ First of all, it is obvious that there is an 
exchange symmetry between Aj and K\ such that 

= (21) 

Secondly, based on the fact that there is a symmetry 
property for the product of U J kz and lf k , z 

( 22 ) 

that is valid no matter which sub-blocks (E + , E“, 0 + , 
and O - ) the wave functions of H 2 0 states belong. With 
this symmetry property, one can show that 

= <-*,-*,■ ( 23 ) 

As a result, we know that some of K , with 
different K\ and K\ are identical. For example, 

^222 - ^ 2 - 2 - 2 . #220 - ^202 = ^ 2-20 = ^ 20 - 2 > aad 
so on. 

On the other hand, the four correlations (2 2 2 0), 
(2—2—2 0), (2 2—2 0) and (2—2 2 0) are identical and the 
other four correlations (2 2 2 2), (2—2—2 2), (2 2—2 2) 
and (2—2 22) are identical. Then, with Equation (11) 
one can conclude that contributions to S 2 , outer, Ac) and 
S 2 , outer, f( r c) from the correlations (2 2 2 0) and 
(2—2—2 0) are the same because not only 
W^2 — ”2-2-2’ ^ Ut a ^ S0 these two correlations are 
identical. Among the eight correlations listed above, 
similar conclusions are also true for the pair of (2 2 2 2) 
and (2-2-2 2), for (2 2-2 0) and (2-2 2 0), and for 
(22—2 2) and (2—2 2 2). There are many similar 


identities valid for other correlations. Based on these 
identity properties, one is able to significantly reduce 
CPU times required when performing numerical calcu- 
lations for the half-widths of H 2 0 lines. 

Besides, with Equation (16) one knows that a 
summation of A n (L\K\K { \ jiZj) over all n equals to 
2L| + 1 if K l =fC l or equals to 0 if K\fK\. For 
examples, the total intensity of the spectrum for Wf 22 is 
5, meanwhile the summation of the Fourier series of 
2 is zero. For the latter, the summations over its 
positive coefficients and over its negative ones are 
equal to 2.5 and —2.5, respectively. This implies that 
when contributions to S 2 ^ outer ,,{r c ) and Stouter, A r c ) 
from the correlations (22-20) are added up, significant 
cancelations would happen between the positive and 
negative coefficients associated with W^_ 2 (t). As a 
result, the net contributions from (22-20) are dramat- 
ically reduced. Numerical results show that in com- 
parison with (2220), the contributions from (22-20) are 
several times smaller. Given the fact that these two 
correlations are completely identical, it is the different 
properties of and Wf_ 22 that play crucial roles 
here in determining which one is dominant. With 
respect to other correlations, similar conclusions are 
also true. Based on these conclusions, one knows which 
correlations are important and which are not. This is 
helpful in our quantitative analyses. 


4.3. Spectra of ^ lKlK {(o;jx) 

In order to exploit the benefits of the ‘Fourier series’ of 
Wf KK (t;./ T ) and to exhibit their intrinsic properties, 
one examines how they vary with members within the 
sets of paired states defined previously in analyzing the 
energy levels and the wave functions. According to 
definitions of the sets, states belonging to the same sets 
are paired with their partners and each of the pairs is 
distinguished by one parameter (i.e., the quantum 
number /). Then, by applying the knowledge on the D 
matrices developed in Section A- 2 of the Appendix, 
one can predict that there are identities between spectra 
associated with two paired states in the same sets, and 
one can also expect that these distribution patterns 
would vary smoothly as the pairs of interest vary. 

We choose J A) as an example. In calculating 

contributions to ReS 2> outer , ,{r c ) or ReS 2 , oute ,,/(r c ), 
these spectra are the most important ones because 
the correlations (1000) and (1002) associated with them 
are the largest among all the correlations. First of all, 
we consider the sets of two paired states and 

jj-n,n+ i with j—jbd, jbd +!,■■■ where n could be one 

among 0,1, We are interested in comparing the 

‘Fourier series’ of IF^/G/n) and ”wo( ( ’J T 2 ) where 
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r i —j — 2n and x 2 —j — 2n— 1. The former is a two 
dimensional array constructed by a>y ri/ v, and 

( 1 00; jx\ ) withy' and x\ (or/', K al , and kf cl ) run- 
ning over all allowed selections where the expression 
for A// r '(100; jx\) is given by 

Ayr’ (100; yn ) = (2/ + 1 )[D(y)-„,^' , , Kt ! 1 «)] 2 . (24) 

Meanwhile, the latter is an array constructed by 
ojyyy and Ajiy^lOO; jxj) withy' and x' 2 over all allowed 
selections where Ay z > ( 1 00; /t 2 ) is given by 

A i% (m-jx 2 ) = (2/ + l)TO_ njI+l /h >t , ; 10)] 2 . (25) 

Because the two original states y}-„, n andy)_„,„ +1 
differ from each other only by k c — n and k c — n + 1 , 
their non-zero coupling partners determined by the 
selection rule must have the same 7' and k' a , but with 
different k' c values differing by 1. For example, as one 
picks a subsidiary state jj,_ n n+l withy' —j+ 1 from the 
allowed coupling list for the original state y}_„,„, one can 
find its corresponding subsidiary one n appearing in 
the allowed coupling list for the original state y}_„, n + 1 ■ 
It is obvious that these two subsidiary states jj,_ n n+x and 
j',_ n n are paired states. Thus, one can conclude that the 
two D matrices D(jj_ nil jj,_ nn+l ; 10) and 
D(jj- n , n +\jj'_ n 10) which appear in Equations (24) 
and (25), respectively, are paired also. Besides, the 
energy difference between the original jj- n ,„ and its 
coupling jj,_ n n+x approximately equals to the difference 
between the originaly)_„ i „ + i and its coupling /)', n . This 
implies that the ‘Fourier series’ of W^^(t\jx) associated 
with y)_„ „ and that associated withy}_„ „ +1 have one 
identical component. By extending this discussion to 
other allowed couplings in their lists, one can conclude 
that these two paired states almost share identical 
components. Similar discussions can be carried out for 
the sets of pairs states j n j- n and y„ +1-/ _„ with 

n = 0, 1, In summary, one should expect that any 

pairs of the states whose j values are above / b d should 
share identity ‘Fourier series’. 

Concerning the analysis of how the spectra would 
vary with pairs within the same sets, one has to rely on 
plotting calculated spectra. For better visualization, 
one can consider the ‘Fourier series’ of W^{t\jx) as 
spectra fV\g 0 (co;jx) consisting of many components. 
Each of their components can be imaged as a 
harmonics with cj a (jx) and A„(100;yr) as its frequency 
and intensity. We would like to note that the terms of 
the spectrum and intensity are used here in describing 
the ‘Fourier series’ of W^“g 0 (t\jx). Readers should not 
confuse these with those in describing resonance 
transitions for the H 2 0 molecule. 

In the following, we select several different sets of 
the pairs and calculate their spectra of W]q 0 (co; jx). 


In Figure 4(a-c), we present the spectra for three sets of 
{jj.oj'yi}, {./'/— 1 , 1 ,7/— 1 . 2)5 and {77-2,2, 77-2,3} with 

j —jhd ■ ■ ■ 2 1 whose k a values are the maximum or 
closer to the maximum and whose wave functions are 
plotted in Figure 2. In the plots, only j values of the 
pairs are explicitly provided. Based on these j values 
and definitions of the symbols in the plots, one can 
easily determine to which of the states the spectra 
belong. As shown in Figure 4(a-c), there are three or 
four branches consisting of main components from 
different states. By comparing the spectra between two 
paired states, one can conclude that their spectra are 
almost identical because two symbols □ and x are 
always overlapped. The higher the j is, the more 
identical their spectra. In addition, by looking at 
structures of the branches one can find that their 
patterns are well organized. This reflects the fact that 
their spectra vary very smoothly as 7 varies. Finally, by 
comparing the plots associated with the different sets, 
one can conclude that different sets have different 
spectral distributions. It is worth reiterating that the 
above claims are valid only for those states whose j 
values are above certain boundaries / b d and the latter 
depend on the sets of interest. 

Similarly, in Figure 5(a-c) we present the spectral 
distributions for another three sets of pairs {joj,jij}, 
{jij-ujij-i}, and { 72, 7-2,73,7-2} whose k c values are 
the maximum or close to it and whose wave functions 
are plotted in Figure 3. We don’t repeat the discussion 
and conclusions here because they are the same as 
those for the previous three sets. In summary, there are 
spectrum pair identities and spectrum pattern similar- 
ities within the sets defined previously. 

There are several others of Wl KK , (®17 t ) with 
different choices of L\, K \ . and K\ . We expect that they 
could differ from each other dramatically. For exam- 
ple, the spectra of W^_ 2 (co;jx) contain components 
with either positive or negative intensities such that its 
total intensity is zero. However, we believe that the 
pair identity and smooth variation features in describ- 
ing W\g 0 (co;jx) within the same sets should remain 
valid for them also because these features result from 
the same sources: the properties of the energy levels 
and the wave functions. 


4 . 4 . The ‘Fourier series’ of W^(t) 

We now present profiles of the ‘Fourier series’ of 
with Z. 2 = 0, 2, . . . which are common for all H 2 0 lines. 
It is obvious that the ‘Fourier series’ of IF < 0 / ’ ) (7) consists 
of only one term with its magnitude equals to 1 and 
a> — 0. Next, we consider the ‘Fourier series’ of \ vf\t) 
which is the most important one among all W^\t) and 
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FREQUENCY (cm -1 ) 



Figure 4. The spectra W^(co;jz) of three sets of the paired states {jj,oJj,i}, and {jj- 2 , 2 , jj- 2 , 3 } with j=jh &, . . . , 21 

whose k a values are the maxima or closer to the maxima. For each of the pairs, spectra of the two paired states are plotted by 
symbols of □ and x in (a) - (c), respectively. Meanwhile, their j values are printed above the symbols □. 


present its spectrum at T— 296 K in Figure 6. As shown 
in the figure, the spectra consists of many components, 
but the one located at <m = 0 is overwhelmingly larger 
than the others. 


4 . 5 . Samples of the ‘Fourier series’ of 

<(') 

We now consider spectra associated with the ‘Fourier 
series’ of K ,(f, j't) We note that although 


are independent of the temperature, 

but W^ h \t) depends on the temperature slightly. As a 
result, the ‘Fourier series’ of 
would vary slightly with T also. 

First of all, it is obvious that the spectra associated 
with t) are the same as that of 

k' ( w ; 7 r i) because the spectrum associated with 
Wq (/) consists of only one component located at a> = 0 
with a unit intensity. With respect to the spectra for 
one expects their main patterns 
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Figure 5. The same as Figure 4 except for three sets of pairs {jojJij}, {jij-ujij-i}, an d {J 2 J- 2 J 3 J- 2 } whose k c values are the 
maxima or closer to the maxima. 


are similar to that of (&>; /r) because the 

dominant component of is located at ®=0 

and its other components are weak. 

In order to demonstrate this claim more clearly, we 
present the spectra of (co; jr) for two paired states 
lOo.io and 10i 10 in Figure 7. As shown by the figure, 
their spectra contain two major components account- 
ing for 90.7% of the total intensity and one minor one. 
Then, we present the spectra of the ‘Fourier series’ 
lF^)(t; jx) f° r them in Figure 8. By comparing 

Figures 7 and 8, it is clear that the whole spectra share 


the main features of W[q 0 (co; jr), but have more small 
structures due to co-adding many weak components 
in wf{oj). 

In the present study, we want to investigate why 
and how the calculated half-widths vary with lines of 
interest. Because the Fourier transforms H l,k,k\l 2 ((ci) 
are common for all lines, the ‘Fourier series’ of 
IF^^i K , (f, jr) contain all the information neces- 

sary to distinguish contributions to ReS^, outer, i O'c) and 
ReS^ outer,/ O'c) for different H 2 0 lines. Given the fact 
that the spectra of wf^(co) are independent of lines. 


Downloaded by [NASA Goddard Institute for Space Studies] at 07:23 23 April 2013 


320 


Q. Ma et al. 




Figure 7. The spectra of w\ a 0 g(a > ; jx) for two paired states 10 o , io and 10i >lo - The former is given by x and the latter by □. 


the only factor to distinguish contributions for indi- 
vidual lines of jf if*—Jh t, are their spectra of 
^L l K t K’ 1 ( com Ji T d and In addition, we 

have shown these spectra represent main features of 
the whole spectra well. Thus, we believe all results and 
conclusions obtained from analysing v (co',jr) 

are similar to those that would be derived if the 


entire spectra were taken into account. Thus, in the 
following analyses we only focus on the spectra of 

Alternatively, one can argue that among all possi- 
ble choices of , many of them are zero. For those 
remaining, one can assume |a^| | because 

the rotational constant of the N 2 molecule is small. 
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Figure 8. The spectra of the ‘Fourier series’ of W^ 0 (t',jr) W^it) at /'= 296 K for two paired states 10 o ,io. and 10i 10 . They are 
given by x and □, respectively. 


Then, by approximating co^jr) (=&> /T/ v + by 
i o a (jz ) {—u>j Z j' z '), Equation (20) can be simplified 

2,outer,i0 r) 

E (2X^+1) 

V L,K, K\ L 2 

x ^A a (L\K\K' l \j i Tj)H LiKiKiLl [(o a (jiT l )\. (26) 

a 

Therefore, we have demonstrated that one can 
use the simpler spectra K , (&>; jz) to analyze 

contributions to ReS 2 , outer, ,{r c ) and ReS^ outerJ (r c ). 
In comparison with Equation (20), the number of 
choices for the a in Equation (26) is much smaller than 
that for the y. 


4.6. The term of S 2 , mid die(r c ) 

With respect to S 2 , m iddie( r c ), it can be expressed as 

B 2 ,middle(j c) ~~ ^ f / [ K] K\ (.//’7(/;T) 

LiK,K\L 2 

/ OO 

dt^(t)F LlKlK[L2 (t) 

-OO 

= E P L \ K\ K\ (jf r fii r i)B L I K, K\ i 2 ( a)- 


In the above expression, P 1a K \ (jf r//)r,) and 
B L\K\K!^L-fj c) are defined by 

p u k, k\ (Jftjjyi) 

- (— 1 ) L ' (2ji + 1 )( 2 j f + IWUdfjijf, 1 L x ) 
x D{jiZijiZu L\K\ )D{j f ZfjfZf\ L\ K \ ), (28) 

and 

PctK^yf/c) — E(2/' 2 + l)(2/2 + 1 )Pji 

iii 2 

x C 2 Uzj 2 P2, 000 )H Li K\K l L 1 (p)j 2 j£) 

— E Pf>(p2)H LlKtK ' iLl (a>f), (29) 

J8 

respectively. In comparison with the expressions for 
Rc.S’ 2 , outer , i(r c ) in Equation (20), the expression for 
S 2 , middled c) i n Equation (27) is simpler because it is a 
product of two factors which are not interwoven. The 
factor of BL,K l K' l L 2 ( r c) is a function of r c and it is 
common for all lines. In contrast, the factor of 
Pr A K, k\ {jf yjiZi) is a constant and its value depends 
on the line of interest. Thus, with respect to a specified 
correlation, how its contributions to S 2t mic jdie( r c) would 
vary as r c varies depends only on and 

meanwhile the magnitude for the line /y f r ; would 
be proportional to the value of Fi A K^K\(jft[jizj). 

As mentioned previously, similar to Re l S’ 2j nuter , ,(r c ) 
and Re S 2 , outer, f( r c)> contributions to S 2 ^ mi ddi e ( r c) from 
individual correlations are also additive. It is not 
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difficult to show there are no contributions to 
S 2 , middled c) from correlations with odd L\ because 
the matrix elements D(jxjx; LK) are zero. In fact, with 
Equation (22) one can show that 

DUrjr- LK) — (~\) r+K D(Jxjx; LK). (30) 

Because we use the II R representation to describe the 
wave functions of H 2 0, the indices K\ and K\ of the 
correlations must be even. As a result. Equation (30) 
means that D( jrjr\ LK) are zero unless L is even. 
Therefore, we know that there are no contributions to 
S 2 , middled <) from all the correlations with odd L\ 
values. 

At this stage we would like to discuss the symme- 
tries associated with P^K^Uf^d and B L{K ^ Ll (r c ). 
First of all, with Equation (22) one can also show that 

D(jxjx\ L — K) = Dijxjx; LK). (31) 

Based on this one can conclude that, as long as 
remains the same, the values of P i A K^K\(jfXfjix) asso- 
ciated with different K\ are identical. The same 
conclusion is also true for K\ . One can conveniently 
express this symmetry as Pi A \K t \\K\\{jfXjjiXi). However, 
there is no exchange symmetry between Aj and K\ in 

P L\K\ K' x (.// TfjjTj) . 

Secondly, we discuss symmetry properties of 
Bl\K\K\l 2 (?c)- With Equation (29), one can conclude 
that symmetries of B^^k^l^c) are determined by 
H l^k i^iohn) and the latter’s are the same as 
those of the correlation function. Therefore, we know 
that B 2 22o( r c) — -#2-2-2o( A-) = B 2 2-2o( r c) = B 2 _ 220 (r c ) 

and B 2222 (r c ) = B 2 _ 2 _ 22 (r L ) — B 22 _ 22 (r c ) = B 2 _ 222 (r e ). 

Similarly, one can express this symmetry as 
^£,|a:i||a:' i |z. 2 (a)- In addition, there is an exchange 
symmetry between K x and K { in B L{KiK ^ Ll {r c ). 

Finally, it is worth mentioning here that with 
respect to making contributions to S 2 _ mi ddi e (r c ), corre- 
lations with K\ ^ K\ and K\ = K\ are comparable. This 
contrasts with S 2 , outer, i{r c ) and S 2 , outerJ (r c ) where 
correlations with K) ^ K\ are much less important 
than those with K\ — K\. As a result, in comparison 
with S 2 , outer,, {r c ) and S 2 , outer j{r c ), how important 
S 2 , middled <!) is would depends on whether K\ — K x or 
K\ ^ K\ for the correlation to be considered. In 
general, S 2 , ou ter, ,{r ( ) and S 2 . outer, /(r,) are larger than 
S 2 , middled c ), but this may not be true for correlations 
with Ki^K\. 

4.7. Simplification in evaluating S 2 , m iMie( r c) 

It turns out that among all B Ll K l K' l L 2 ( r c), the most 
important ones are B 222 o(r c ) and B 2222 (r c ). Then, one 


can adopt a simpler approximate formula to calculate 

B 2 , middled c) 

S 2 ,middle( r c) — A >'jLn P 222 (./) T.// Fd [ ^2220 fr c) + B 22 2 2 ( r c )]> 

(32) 

where the factor of 4 appearing on the right side of 
Equation (32) results from the four identical contribu- 
tions for each B 222 o( r c) and .82222(A)- From Equation 
(32) one can draw two important conclusions about 
S 2 , middled c)- First, the pattern of S 2 , middled c) (i-e., as a 
function of r c as it varies) is mainly determined by 
B 22 2 xky c) + B 2222 {r c ). Second, for specified H 2 0 lines 
their contributions to S 2 ^ mic idie( r c) are roughly propor- 
tional to values of F 222 {.if x Ji r ,)- Thus, with these 
values, one is able to easily judge which lines would 
have large magnitudes of S 2t mit jdte( r c) and which lines 
would have small ones. Therefore, values of 
F 222 (y/T//,-r,) play a crucial role to distinguish 
S 2 . middled c ) for different lines. 

Based on the potential model used in updating 
HITRAN 2008 [17], the 20th order cut-off and the 
‘exact’ trajectory model, we have derived the function 
of 4-v/27r[R 2 22o(A) + -62222(A)] at T = 220, 296, and 
340 K and present results in Figure 9. As shown in 
the figure, values of 4v / 27r[5 2220 (r f .) + B 2222 {r c )\ 
decrease as r c increases from the starting point A.min- 
Meanwhile, by comparing its magnitudes at different 
temperatures, one can conclude that the magnitudes 
slightly decrease as T increases. 

On the other hand, one can easily evaluate values of 
P 222 {jf r tii a) for all H 2 0 lines of the pure rotational 
band listed in HITRAN from Equation (28). In 
contrast with 4V2n[B 222 o(r c ) + 62222(A)], ^(i/f^A) 
is a constant, it does not depend on the temperature, 
and its values are independent of the potential and 
trajectory models used in calculations. When values of 
P 2 22 Uf r t)i r d for all H 2 0 lines are available, it becomes 
very easy to obtain the corresponding S 2 ^ midd i e {r c ) 
terms approximately for individual lines of interest 
because one only needs to multiply 
4V2 tt[ 6 2 22 o(a) + 62222(A)] by their P 2 2 2 {jfX l j i x 'i) values. 

In order to exhibit intrinsic properties of 
P 222 Uf T jii r i), one needs to calculate values for lines 
belonging to the individually defined groups. 
In Section 5, we will provide their definitions in 
detail. At present, we only note that each of the groups 
consists of a set of paired lines in a specified branch. 
As an example, we choose three groups, {jj qj. 4-jnj,, 
jftjf jiOji h {jfijf-i ji 2 j l -\,jf 2 jf-\ 1}, and 

{jfljj -3 *-ji 2 j t -\,jf 4 jj ~3 in the R branch of 

the pure rotational band and present their calculated 
values in Figure 10. As shown in the figure, for two 
paired lines whose j values are not less than either of 
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Figure 9. A plot to show profiles of 4\/2jr [5 2 22o(/c) + Bi 2 ii( r c)\ at 7’=220, 296, and 340 K. They are plotted by solid, dashed, 
and dash-dotted curves, respectively. 



3 4 

5 6 7 8 9 

10 

11 12 13 14 15 16 17 

J i 

18 19 

20 21 

Figure 10. Calculated values 

of R 222 O/W 1 ) 

for 

three groups of FFO 

lines 

Ujoj, * jnj,-jn,i, ' M/. 


{jfijf-i jajt-iJfijf-i i}, and {jfsjj-i ■*- jaj,-i Jftjj-i jnj ,- il in the R branch listed in the pure rotational band 

of HITRAN. They are plotted by three pairs of symbols {x, A), {+, Oh and {*, □}, respectively. 


the /bd associated with their initial and final states, their 
PmUfyiPd values are always identical. More specifi- 
cally, for these three groups, as long as their members’ 
ji values are above 7, 10, and 15, respectively, 
their paired lines have identical Pz22(Jf x fji r d values. 


Meanwhile, by looking at how their values vary as 
the pair of interest varies within the same groups, 
one can conclude that the variations with /,■ 
are very smooth. Finally, by comparing values 
of P 222(7/ Td A/) th e same ji among different 
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Figure 11. The same as Figure 10 except for three groups {j f2 j f ~i < jnj r h\jt i * jaj, \J/y f 2 and 

UfiJr-i jajt-iJftjj-i ■<-jaj i - 2 } in the Q branch. 



groups, one can conclude the values could change 
dramatically. 

As another example, we present results of 
PiiiUfyii*,) for three groups ^jnj„ 

jfijj-i {jfijf-2 Jfijf-2 1} 5 an d 

{7/3 ^-3 jfijf-i ■*— jaji-i} in the Q branch in 

Figure 1 1 . We don’t repeat the discussion here because 
the main conclusions are the same as above. We only 
note that for paired lines in these three groups, the pair 
identity of PmUfyi^i) is valid starting from /, = 10, 
13, and 16, respectively. 


5. Categorizations of H 2 0 lines and profiles of 
ReiS 2 (r t .) for lines within the same groups 

We have shown that by dividing FI2O states into the 
sets of paired states, there are the pair identity and 
smooth variation properties applicable for the energy 
levels, the wave functions, and the spectra 
”LiK I K'( C0 ’j T ) °f states belonging to the same sets. 
This knowledge provides a useful hint for how to 
divide the H 2 0 lines into groups such that for lines in 
the same groups their calculated N 2 broadened half- 
widths have similar properties. It is obvious that in 
specifying the lines, one has to categorize their initial 
and final states simultaneously. First of all, in order to 
reduce one variable from /, and jf, one divides lines into 
the P, Q, and R branches. Then, by determining which 
of the sets of paired states their initial and final states 


belong to, one categorizes those lines in the same 
branches. 

For example, a group in the R branch {jf— y) + 1) 
consists of a set of pairs of lines {jjoj f •*- 7; 17,, 
with ji =0,1, 2,.... In this group, two 
lines with the same j t are called paired lines because 
their final sates jj^j, and jf\j f are paired partners in the 
same set and their initial states juj, and 7/07, are also 
paired partners in the same set. In this group, the pairs 
of lines are determined by only one variable j, and the 
interesting pairs are those with /,■ above the boundary 
/ bd (=7). Due to similar properties shared by their 
energy levels and wave functions, one expects that their 
calculated N 2 -broadened half-widths must somehow 
follow a similarity too. In fact, for each of the paired 
lines in this group we have explicitly shown the spectra 
W\ g 0 (co;jT) associated with their initial and final states 
in Figure 5(a) and provided their values of P222Uf r jft T d 
in Figure 10. Given the fact that the former mainly 
determine their S 2 , outer, and S 2 , outer , / (a) terms and 

the latter determine their Sf m iddie( r c), one can be pretty 
sure that their half-widths must have similar proper- 
ties. More specifically, as long as 7, > 7, the paired lines 
in this group should have the same half-width values. 
In addition, these values should vary smoothly with 
the pairs. 

In order to demonstrate this claim, we present 
calculated Rc.S' 2 (a) at T — 296 K for lines with 
ji—1, 8 , . . . , 19 in this group in Figure 12. The calcu- 
lations are carried out with the modified RB formalism 
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3.5 3.55 3.6 3.65 3.7 3.75 3.8 3.85 3.9 3.95 4 


T (A) 

Figure 12. Calculated Rc.S’ 2 (r c ) at T= 296 K for lines with /,= 7, 8, ... , 19 in the group of jpj. 4- j n j, and jf\j, 4— j i 0 j. in the R 
branch. For each of the paired lines, their ReS^Ov) are plotted by a solid and a dotted curves. Results associated with different j t 
values are plotted by different colors and lie from top to bottom in the figure. 


and based on the potential model used in updating 
HITRAN 2008 [17]. In addition, 132 correlations have 
been taken into account and the ‘exact’ trajectory 
model was used. As shown in the figure, for each of the 
pairs their ReS^iy) are always identical and the 
profiles of Rc.SA/y) associated with different /,■ values 
change very smoothly. Based on this, there are pair 
identity and smooth variation properties for their 
calculated half-widths because the latter are completely 
determined by these functions of ReS 2 (r £ ). We will 
return to this subject latter. 

Similarly, one can define other groups in the R 
blanch such as { jfnjf—n * jin- 1 ./, nijfn- 1 ./, n * jinjj—n} 

with n— 1,2, In defining these groups, values of k c 

for their initial and final states are closer to the 
maxima. Meanwhile, we can define groups in the R 
branch where values of k a are equal or close to the 
maxima, such as groups consisting of paired 

lines { jfjj—n,n 4 jij)—n,n+] ijf/j—n,n+] 4 jiji—n,n\ With 

n = 0,l, As long as the selection rule is allowed, 

more groups based on many other combinations of 
paired lines are accessible. In addition, there are groups 
consisting of un-paired lines. For example, 
jfjf- 1,2 4 — jioji with j)— odd are considered as un- 
paired lines because one cannot find other lines to 
pair both their final and initial states. 

For the P and Q braches, one can perform similar 
procedures. As another example, we can define a group 
in the Q branch consisting of a set of paired lines 
jf2jf-\ 4 - jnji and jfijj-] with = 1 , 2 ,.... 


We present calculated Re5 , 2 (7' c ) at T= 296 K for pairs 
of lines with ji — 9, 10, ... ,20 in this group in Figure 13. 
As shown in the figure, profiles of these functions of 
ReS 2 (r c ) have similar features as the group in the R 
branch given above. 


6. Pair identity and smooth variation rules for the 
half-widths 

After categorizing H 2 0 lines based on properties of 
energy levels and wave functions of their initial and 
final states, we are ready to establish rules applicable 
for the calculated half-widths for lines within individ- 
ual groups. 

We select six groups of lines in the pure rotational 
band and present their calculated N 2 -broadened 
half-widths in Figure 14(a f). They are three groups 

Ufijf Aoy,}j jaji-\ , 

if and {jf2jf-2 ^jnjiJfijf-2 «- jmj,} in 

the R branch, two groups of 

}f 2 j r \ and {y/ 2^-2 } 

in the Q branch, and one group of 

{jf 2 jf -2 jiuJfSjf -2 jioji } in the P branch. As 
shown in the figures, paired lines with their j, values 
above boundaries in all six groups have almost 
identical half-widths. In addition, by comparing dif- 
ferent pairs in the same groups, their half-widths vary 
smoothly as their /, values vary. In general, the 
boundaries of the groups of lines are determined by 
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(A) 

Figure 13. The same as Figure 12 except for lines with /,-= 9, 10, . . . , 20 in a group of jf2j, \ ■*- jn j, andjfij i <— j^j. in the Q 
branch. 


states with the smallest of \k a — k c \. With Table 2, one 
can predict these boundaries and they are 7, 10, 13, 10, 
13, and 14, respectively. We note that for lines in the 
P branch, their jf are less than /,. As a result, the 
boundary requirement applies for jf (i.e. jf> 13) and 
the latter means /,•> 14. 

Thus, the pair identity and smooth variation rules 
applicable for calculated half-widths for lines involving 
high j states within the groups can be established. The 
higher the j b the firmer these two rules hold. 
As explained above, the very origin of these two rules 
is the properties of the energy levels and the wave 
functions of H 2 0 states. No matter what kinds of 
potential and trajectory models are used in calcula- 
tions, theoretically calculated N 2 -broadened half- 
widths must follow these rules. 

Furthermore, by scrutinizing the whole process pre- 
sented above, one cannot find any reason why 
the analyzing process and the final conclusions 
described above for the widths are not applicable for 
calculated pressure-induced shifts. In fact, with the 
modified RB formalism, the expression for the shift is 
given by 

/H-oo 

vf(y)dv I 2nbdb sin((S\(b)) 

Jo 

+ Im(5 2 (h)))e- Re<S2(i}> . (33) 

For lines in the pure rotational band, by replacing 
the integration over the velocity by the averaged 


s M RB = r 

2jtc J 0 


velocity the above expression can be simplified as 



sin[ImS 2 (r f )]e ReS2<rt W c . 


(34) 


Because Sf midd / e (r c ) is real, Im5 2 (r f ) consists 
of only two components -Im S 2 , outer, iVc) and 
Irn S 2 , outer, f(r c ) where the minus sign of ImS 2 , outer, ,{r c ) 
results from the fact that S 2 (r c ) = S 2 , outer, t {r c )* + 
*^ 2 , outer, f (r c ) + S 2 , middled c ) [ 11 ]- The expression for 
Im S 2 , outer, i( r c) is the same as Equation (9) for 
ReS 2j outer, /(a) except for a replacement of 

H L\K\ k’ x Li (^a) by the Cauchy principal integrations 
Iuk^iAco) defined by 

1 f +co , 1 

luKiK'.Ltiw) — P / da)'— kf [ lKlK 'L 2 (oJ), 

71 J—OO &> — CO 

(35) 


where P means the principal part. Similar to the 
function of Il x k,k\l 2 {u>) are common for 

all lines. As a result, all the discussions given above are 
applicable for the calculated shifts and the same two 
rules must be valid also. 

One can extend the discussions to other tempera- 
tures for the H 2 0-N 2 system. One can carry out similar 
analysis for other systems involving different pertur- 
bers such as H 2 0-0 2 and H 2 0-H 2 0. It seems that 
there are no any essential differences occurring among 
these different cases. Therefore, the two rules are not 
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Figure 14. Calculated N 2 -broadened half-widths for paired H 2 0 lilies in six groups defined in the pure rotational band. Three 
groups in the R branch are plotted in (a)-(c), two groups in the Q branch are in (d)-(e), and one group in the P branch is in (f). 


only applicable for the N 2 -broadened half-widths, but 
also for the air- and self-broadened half-widths and 
their corresponding temperature exponents. 


7. Discussion and conclusions 

Based on this understanding acquired in the present 
study, one can develop useful ways to improve the 
accuracy of the spectroscopic parameters in databases. 


First of all, by applying the identity rule to two paired 
lines, one is able to pick out those values listed in the 
database that differ appreciably from one another. 
One can be fairly sure that those data picked out con- 
tain large errors or, for instance, came from different 
theoretical or experimental sources. Meanwhile, by 
applying the smooth variation rule for lines in the same 
groups, their individual errors could be partially aver- 
aged out using proper smoothing procedures. Thus, 
one is able to obtain better values with less uncertainty. 
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For people who are in charge of databases, the 
present work provides a helpful guide and a useful tool 
in their developing practices. It is wise to consider that 
all lines in each of the transition bands belong to a well 
organized network. Acting as one member of the 
network, the half-width value of an individual line is 
not self-contained because it has to be somehow 
connected with other values associated with those 
lines, especially its paired line, who belong to the same 
group as the line of interest does. The same statement 
is also true for other spectroscopic parameters. Given 
the fact that data listed in the databases may come 
from several independent sources, to take a balance 
procedure becomes even necessary. Meanwhile, by 
applying the two rules to each of the groups of lines 
thoroughly, one can screen the database and identify 
lines whose half-width values (or other spectroscopic 
parameters) that behave strangely. These outliers 
should be carefully checked, or one could consider 
the option to replace them by smoothed values 
obtained from their neighboring lines. After removing 
the errors, they are able to improve the accuracy of the 
database significantly. 

Experimentalists can also use the pair identity and 
smooth variation rules to check their measured spec- 
troscopic data. If their results scatter wildly, it 
definitely indicates contaminations or errors existing 
in their data. In addition, the smooth variation rule 
provides a sound support and is a useful guide for their 
developing data smoothing procedures. Again, one has 
to keep in mind that although the smoothing practice 
for data with high j states in the same groups is well 
justified, one cannot apply it to those lines with low j 
states, otherwise the procedures could introduce arti- 
ficial errors. Finally, with respect to the determination 
of measurement priorities, it is prudent to choose 
several lines with high js in the same group because 
several good measurements are enough to determine 
others in the same group. This implies that it is not 
necessary to measure all lines in the same group for the 
completeness purpose. 

In summary, it is the properties of energy levels and 
wave functions of the H 2 0 states that play an essential 
role in establishing the pair identity and the smooth 
variation rules and these rules reflect natural correla- 
tions between the inputs and the outputs for a whole 
system consisting of an absorber H 2 0 molecule 
immersed in bath molecules and radiation fields. 
Flowever, due to the complexity of dynamical processes 
happening inside the system, these correlations are very 
difficult to be fully grasped. This implies that it is 
unrealistic to set as a goal to find common rules with 
which the outputs can be well monitored from the 
inputs unless one narrows the variation range of the 


inputs. Thus, it is the categorization of the H 2 0 lines 
that enable one to narrow the variation ranges of the 
inputs first. In fact, after completing the categorization 
procedures for each of the P, Q , and R branches, there is 
only one independent variable (i.e. the initial quantum 
number /,) left to distinguish lines of interest. Then to 
establish the rules valid within the individual groups 
becomes possible. Therefore, these two rules have two 
characteristics: they are natural, but they are local. 
Local here means the rules work for each individual 
groups and are valid only for its members above certain 
boundaries. Finally, while the present study was carried 
out in the pure rotational band of H 2 0, one can extend 
the study to other H 2 0 bands and establish similar rules 
there also. Of course, one has to analyse the properties 
of the energy levels and wave functions of the states 
involved in these transitions first. 


Acknowledgements 

Two of the authors (Q. Ma and R.H. Tipping) acknowledge 
financial support from NASA under grants NNG06GB23G, 
NNX09AB62G, and FCCS-547. Q. Ma wishes to acknowl- 
edge financial support from the Biological and 
Environmental Research Program (BER), US Department 
of Energy, Interagency Agreement No. DE-AI02-93ER61744 
and financial support from NASA under grant 
NNH08ZDA001N-ACLAB. This research used resources 
of the National Energy Research Scientific Computing 
Center, which is supported by the Office of Science of the 
U.S. Department of Energy under Contract No. DE-AC02- 
05CH1 1231. 


References 

[1] L.S. Rothman, I.E. Gordon, A. Barbe, D.C. Benner, 
P.F. Bernath, M. Birk, V. Boudon, L.R. Brown, 
A. Campargue, J.-P. Champion, K. Chance, 

L. H. Coudert, V. Dana, V.M. Devi, S. Fally, 

J.-M. Flaud, R.R. Gamache, A. Goldman, 

D. Jacquemart, I. Kleiner, N. Lacome, W.J. Lafferty, 
J.-Y. Mandin, S.T. Massie, S.N. Mikhailenko, 

C. E. Miller, N. Moazzen-Ahmadi, O. Naumenko, 
A.V. Nikitin, J. Orphal, V.I. Perevalov, A. Perrin, 
A. Predoi-Cross, C.P. Rinsland, M. Rotger, 

M. Simeckova, M.A.H. Smith, K. Sung, 

S.A. Tashkun, J. Tennyson, R.A. Toth, A.C. Vandaele 
and J. Vander Auwera, J. Quant. Spectrosc. Radiat. 
Transfer 110 , 533 (2009). 

[2] I.E. Gordon, L.S. Rothman, R.R. Gamache, 

D. Jacquemart, C. Boone, P.F. Bernath, 

M.W. Shephard, J.S. Delamore and S.A. Clough, 
J. Quant. Spectrosc. Radiat. Transfer 108 . 389 (2007). 

[3] D. Robert and J. Bonamy, J. Phys. (France) 40, 923 
(1979). 


Downloaded by [NASA Goddard Institute for Space Studies] at 07:23 23 April 2013 


Molecular Physics 


329 


[4] Q. Ma, R.H. Tipping and R.R. Gamache, Mol. Phys. 

108 , 2225 (2010). 

[5] Q. Ma, R.H. Tipping and N.N. Lavrentieva, Mol. Phys. 

109 , 1925 (2011). 

[6] G.W. King, R.M. Hainer and P.C. Cross, J. Chem. 
Phys. 11 , 27 (1943). 

[7] H.C. Allen and P.C. Cross, Molecular Vih-Rotors 
(Wiley, New York, 1963). 

[8] W. Gordy and R.L. Cook, Microwave Molecular 
Spectra (Wiley, New York, 1970). 

[9] R.N. Zare, Angular Momentum (Wiley, New York, 
1988). 

[10] R.J. Barber, J. Tennyson, G.H. Harris and 
R.N. Tolchenov, Mon. Not. R. Astron. Soc. 368, 1087 
(2006). 

[11] R. Lynch, R.R. Gamache and S.P. Neshyba, J. Chem. 
Phys. 105 , 5711 (1996). 

[12] R.R. Gamache and J. Fischer, JQSRT 78, 289 (2003); 
78. 305 (2003). 

[13] R.R. Gamache and J.-M. Hartmann, J. Can. Chem. 82 , 
1013 (2004). 

[14] R.R. Gamache, J. Mol. Spec. 229 , 9 (2005). 

[15] Q. Ma, R.H. Tipping and C. Boulet, J. Quant. 
Spectrosc. Radiat. Transfer 103 , 588 (2007). 

[16] Q. Ma, R.H. Tipping and C. Boulet , J. Chem. Phys. 
124 , 014109 (2006). 

[17] R.R. Gamache and A. Laraia, J. Mol. Spectrosc. 257, 
116 (2009). 


Appendix 

A. I. Symmetry assignments of H 2 0 wave functions 
for paired states 

We consider the symmetry assignments for the sets of paired 
states and 7/-«,«+i where n = 0, 1, ... whose wave 

functions in the 1 R representation are given in Figure 2. 
It is obvious that the paired states have the same k a values, 
but their k c values differ from each other by 1 . This implies 
that they have identical evenness or oddness for k a , but have 
opposite evenness or oddness for k c . Thus according to 
Table 1, no matter whether their j values are even or odd, 
their symmetry assignments by the superscripts + and — in 
the I R representation must be opposite. In terms of their 
coefficients Uf, the coefficients of one state must be an even 
function of k and the coefficients of its paired partner must 
be an odd function of k. In Table A-l, we list the assignments 
for these pairs. 


Table A-l Symmetry assignments for pairs of and 

jj-„,„+ i in the I R representation. 


jjfii jj, 1 jj- 1,1 j jj- 1 ,2 jj- 2,2, jj- 2,3 

y' = even E + , E~ 0 + , CT E + , E“ 

y = odd 0 + , cr E+, E" 0 + , CP 


Similarly, we consider the symmetry assignments for the 
sets of paired states jnj-n and j n+ ij-„ where 
n = 0, 1, . . . whose wave functions in the III R representation 
are given in Figure 3. In this case, the paired states have the 
same k c and their k a differ from each other by 1 . Then, one 
can conclude that with respect to their wave functions 
derived in the III R representation, one must be an even 
function assigned by the superscript + and its partner must 
be an odd function assigned by the superscript—. In terms of 
their coefficients U' kr , the conclusion is the same as that for 
the sets of and 1 drawn above. In Table A-2, we 

list the assignments for them. These symmetry properties are 
important and we will be back to this subject later. 


A. 2. Important quantities in determining couplings 
between H 2 0 states 

As shown in Section 4, in the RB formalism there are two 
quantities which play crucial roles in determining how the 
Lorentzian half-widths would vary with lines of interest. For 
two H 2 0 states jz and j'z' which can be coupled by the 
irreducible tensor L and its subsidiary index K, the quantities 
cojTj'T' defined by 

cojrfe = [£ <H2 ° , (7t) - £ <H2 ° ) (7 V)]/* (A-l) 

represent their energy differences and the quantities 
D(jzj'z'\ LK) called as the D matrices represent their 
couplings; and this is defined by 

DUrj'z '; LK) = ]T(-1 ) ku iX-Ke C(jf L. kK - kK). 

k 

(A-2) 

As shown by Equation (A-l) and (A-2), with the energy 
levels and the wave functions of these two states, one can 
easily calculate the energy differences and the D matrices. 

However, it turns out that exploiting the properties of 
these two quantities is fruitful in answering a fundamental 
question why and how the calculated N 2 -broadened half- 
widths and pressure induced shifts would vary with lines of 
interest. Based on our analysis of the properties of the energy 
levels and the wave functions of the H 2 0 described 
previously, one can find the properties of these two quantities 
accordingly. 

Let us consider two pairs of the states {j z\,j z 2 ) and {j' 
z\,j' z' 2 ). We assume j Z\ and / z 2 are paired and j' z\ and/' z' 2 
are paired. By picking one member from each of the pairs, 
one gets two combinations, for example, a combination of j 
zi and/' z\ and a combination of j z 2 and/' z 2 . We assume 
that the two states in the same combinations can be coupled 
by L and K where L is the irreducible tensor rank and K is its 
associated index [4,16], and we would like to compare the 
energy differences and the D matrices between these two 


Table A-2. Symmetry assignments for pairs of j n j-„ and 
j„ + \j-n in the III R representation. 


JO.p J 1 J j\ J- 1 5 Jl,j - 1 J2,j—2i Jl.j-2 

7= even E + , E~ CT, 0 + E+, E" 

7 = odd Q~, Q + E + , E“ CT, 0 + 
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Figure A.l. A plot to show properties of D(jk a ,kj k ' k\LK) with L= 1 and K= 0. Calculated values of 1 - 0 ( 7 //, , ; 10) | and 
I D(jj t \jj, 0 : 10) | with / =7+1 and 7 = 3, ... 26 are plotted by x and A, respectively. Because values are positive for even 7 and 
become negative for odd 7 . In order to show how their magnitudes vary with 7 more clearly, their absolute values are plotted in 
the figure. Similarly, a set of pairs | D(jjojL\ j ; 10) | and | 0 ( 7 / 17 '', , 10) | with /' = 7 whose values are negative for even 7 and 

positive for odd 7 are plotted by + and <). Finally, another set o f pairs | D(jj,ojj',v 10 ) I and I D Uj,dj>,o'’ 10 ) I with / =j~ 1 whose 
values are negative for even 7 and positive for odd 7 are plotted by * and □. 


combinations. First of all, it is obvious that there are 
identities of their energy differences 


w j +/+ ~ 


(A-3) 


approximately valid within certain accuracy tolerances when 
both 7 and 7 ' are above corresponding boundaries defined for 
their sets. 

On the other hand, for the two O matrices D(jz\j'z [ ; LK) 
and 0(7+//; LK), one has shown that | U J kzi \ & \U kz , | and 
\U' kz ,\ & | U' kT , | . In addition, as shown in Section A.l, the 
evenness or 2 oddness of U J kTi and U kr n over k must be 
opposite and this is true for lJ kz , and \J kx , also. Thus, it 
becomes clear that 


Tp T P R: Tp TP 
U kri U k-Kzi U kz 2 U k-Kx^ • 

As a result, one is able analytically to verify 


(A-4) 


0 ( 7 + 7 / ; LK) « Dijzjj'p- LK). (A-5) 


Thus, one can draw an important conclusion that the 
coupling by L and K between two states approximately 
equals the coupling between their paired partners. The higher 
the / is, the more equal these couplings become. We call these 
two D matrices as the paired D matrices later. Besides, there 
is another interesting feature of the paired D matrices. This 
feature represents a behavior of their values as 7 varies. 
In order to exhibit it, one has to find their values by 
performing numerical calculations explicitly. 

It is worth mentioning that Equation (A-4) is applicable 
for the coefficients derived accordingly in the I R or III R 
representation, but not in the II R representation. This 
implies that we have verified Equation (A-5) analytically in 


the 1 R and III R representations, but not in the II R 
representation. On the other hand, in developing the line 
shape formalism with the correlation functions described in 
Section 3, we prefer to choose the II R representation in 
which the symmetry axis of the H 2 0 molecule lies along the z 
axis of the molecular fixed frame. With this choice, one is 
able to exploit more symmetries of H 2 0 to reduce a number 
of the correlation functions required to be evaluated 
significantly. Because Equation (A-4) is not valid, to verify 
Equation (A-5) analytically in the II R representation is 
not as easy as shown above. In the present study, we do not 
pursue the analytical verification, rather we carry out 
numerical verification. In addition, numerical results are 
needed to exhibit another feature of the paired O matrices. 

For later convenience, we rewrite D( jzj'z'\ LK) as 
D Uk a ,kj' K , k p LK where k a -k c = z and kf a - kf c = zf . 
We note “that there are rules to determine whether there are 
non-zero couplings between two states j ka kc and j' k , k , or not. 
For example, in the II R representation the rule “associated 
with L = 1 and K = 0 is similar to the selection rule of the 
H 2 0 transitions; Ay (=/'— j) = 0, ± 1, Ak a (=k' a — k a ) = 
±1, ±3, . . . , and A k c (= k' c — k c ) = ±1, ±3, .... The rule for 
L = 2 and K=0 is A/=0,±l,±2, Ak a = 0, ±2, ±4, . . . , and 

A k c = 0, ±2, ±4, For other L and K values, the rules 

become more complicated. 

As an example, we consider a set of two D matrices 
{ L>(Jjfljj, | ; \0), D(j jA j', 0 -, 10)) with/ =y‘+l. It is obvious 
that these two are paired D matrices because j j 0 of the first D 
and 7 / of the second D are paired states, and / y , in the fist D 
and 7/0 in the second D are paired. In addition, there are 
non-zero couplings both between j j0 and 7 / 1 of the first D 
and between j j X and 7 y 0 of the second D. We have calculated 
these two paired D matrices for values of for 7 = 3, ... 26 and 
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Figure A. 2. The same as Figure A. 1 except for L = 2, K= 2, three sets of the paired D matrices ; 22), ; 22) J with 

j'=j+ 1, v 22), ; 22) j withy' =y, and {D(j 0 J[j,; 22 22)) withy' =y- 1. 


found their values are almost identical. However, it turns out 
that the values are positive for / = even and become negative 
for y = odd. In order to show how their magnitudes vary 
with j more clearly, we plot their calculated magnitudes in 
Figure A-l. Similarly, we select another two sets 
- 10),-D(y),iy‘v_i_ 2 ; !0)) with j'=j and 
{DUjfiJyX’ , 0 ; 10)1 with j' =j— 1. Both of them 

have non-zero couplings and both of them consist of two 
paired D matrices. Again, because their calculated values 
change signs for even j and odd j alternatively, we plot their 
magnitudes in Figure A.l. Based on these results, we can 
draw two conclusions. First of all, values of the paired D 
matrices are almost identical. This implies that we have 
verified Equation (A-5) numerically for these sets. 
In addition, we would like to note that for these three sets, 
the identity of paired D matrices becomes valid starting from 
j— 3, 5, and 4, respectively. These j values are consistent with 
those y bd values listed in Table 2. Readers may wonder for the 
third set, why the identity starts from y = 4, but not /= 3. It is 
easy to explain this because for this set, j' =j — 1 and the 
smaller of j and j ' must be not less y bd . Secondly, as shown in 
the figure, one can conclude that for each of these sets, their 
magnitudes vary very smoothly as y varies. Variation of the 


magnitudes with j smoothly is just a feature of the paired 
D matrices we are looking for. 

As another example, instead of choosing states asso- 
ciated with large k a values we consider sets of the paired 
D matrices whose states are associated with large k c values. 
In addition, with respect to L and K we consider a case of 
L = 2 and K= 2. More specifically, we select 22), 

D (pjJoj'’ 22)1 with j'=j+ 1, 22), 

DUijjy’-i '■ 22)} with j'=j, 22), D{j\jjoj,; 22)} 

with j' —j — 1. We present their calculated magnitudes with 
y = 5, ...26 in Figure A. 2. As shown in the figure, the two 
features exhibited previously remain the same. We do not 
repeat the discussion here, other than to mention that in 
comparison with the previous three sets, the identity of two 
paired D becomes valid starting from a higher j value. For 
these three sets, it starts from j =7, 10, and 8, respectively. 
These starting values are consistent with those y bd values 
listed in Table 2. 

With respect to other sets, we believe the same conclu- 
sions mentioned above are applicable there also. As shown in 
Section 3 these two features of the paired D matrices are very 
helpful in analyzing values of calculated half-width for the 
H 2 0 lines. 


